Quasi-two-dimensional non-collinear magnetism in the Mott insulator Sr2F2Fe20S2 
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The magnetism of Sr2F2Fe20S2 was examined through neutron powder diffraction and thermo- 
dynamic and transport measurements. Quasi-two-dimensional magnetic order develops below Tjv = 
106(2) K with an in-plane correlation length exceeding 310 A and an out-of-plane correlation length 
of just 17(3) A. The data are consistent with a two-k structure with ki = (1/2, 0, 1/2) and k2 = (0, 
1/2, 1/2) and an ordered moment of 3.3(1) oriented along the in-plane components of k. This 
structure is composed of orthogonal AFM chains intersecting at super-exchange mediating O sites. 
Density Functional Theory (DFT) also points to this structure and a narrower Fe 3d band than for 
the iron pnictides from which electronic correlations produce a Mott insulator. 



PACS numbers: 75.25.-j, 75.30.-m, 71.20.-b 

Quasi-two-dimensional magnetic materials near the 
metal-insulator transition have been the focus of in- 
tense research since the discovery of high Tf, supercon- 
ductivity in layered copper oxides. The discovery of 
superconductivity in the chemically distinct iron pnic- 
tides has brought attention to the potential for novel 
electronic phases in layered iron-based compounds at 
the local-to-itinerant moment boundary. Here we show 
that Sr2F2Fe2 0S2 is a Mott insulator with quasi-two- 
dimensional non-collinear magnetic order. The unusual 
magnetic structure features two wave-vectors ki = (1/2, 
0, 1/2) and k2 = (0, 1/2, 1/2) associated with perpen- 
dicular AFM spin chains that share super-exchange me- 
diating oxygen sites. 

Sr2F2Fe20S2 is structurally related to La203Fe2X2 (X 
= S, Se, "2322") which was first reported by Mayer et 
apThis is a a layered, tetragonal structure consisting of 
stacked [La202]^^ and [Fe20X2]^~ sheets, composed of 
edge-sharing La40 tetrahedra and face-sharing Fe02X4 
octahedra, respectively. The Fe2 0X2 sublattice forms a 
checkerboard spin lattice, with O and X mediated Fe- 
Fe super-exchange interactions. Insulating behavior and 
long range antiferromagnetic (AFM) ordering at 105 K 
(X = S) or 93 K (X = Se) was reported in these sys- 
tems. Density Functional Theory (DFT) indicates these 
compounds are Mott insulators.l^ 

Very recently, the isostructural compounds 
La2 03C02Se2,'^ R203Mn2Se2 (R = La, Ce, Pr),EHIl 
R203Fe2S2 (R = Ce, Pr)^ and R203Fe2Se2 (R = 
La-Sm)P^ were reported. All "2322" materials known so 
far are insulators or semiconductors with AFM order.^J 
While no band structure calculations were reported 
for the Mn compounds,'^'^ transport and magnetic 
properties similar to those of La203Fe2S2 indicate a 
Mott-localization picture may be applicable. 

By replacing [La202]^^ sheets with different 2+ va- 



lence spacer layers, several other related compounds were 
synthesized, including A2T2OX2 with A = SrF, BaF 
("22212"_p] or Na ("2212"),E3 and T = Mn, Fe and X 
= S, Se.12l For a view into electronic correlations in this 
family of materials, we examine the magnetic structure 
of Sr2F2Fe20S2 finding it distinct from anything seen in 
the copper oxides or the iron pnictides. 

Poly crystalline Sr2F2Fe20S2 was synthesized by solid 
state reaction, using SrF2, SrO, Fe and S powders as 
starting materials. The mixture was pelletized and an- 
nealed in vacuum at temperatures between 800 °C and 
850 °C. Neutron powder diffraction data were collected 
for temperatures between 12 K and 300 K on the time- 
of-flight POWGEN powder diffractometer at the Spalla- 
tion Neutron Source (SNS) at Oak Ridge National Lab- 
oratory. In order to focus on the chemical and mag- 
netic structures respectively, the data were acquired us- 
ing two different instrumental configurations with wave- 
length bands of 1.066 A-2.132 A (Fig. [Tji) and 4.264 
A -5.33 A (Fig. [2]). The room temperature diffraction 
pattern (red points. Fig. [l^) was refined to a tetrago- 
nal structure (Fig. [Tja) using fullprof.^^ The lattice 
parameters a = 4.0362(5) A and c = 17.9915(1) A are 
consistent with the previously reported values.^^ The re- 
fined atomic positions and isotropic thermal factors are 
provided as supplemental material.'^ Minute amounts of 
impurity phases (0.84% FeF2 and 0.39% FeF3) were also 
detected. 

For temperature below Tm — 106(2) K, additional 
peaks (green tick marks. Fig. [2j3 and c) appear in the 
neutron diffraction pattern. These correspond to a prop- 
agation vector k ~ (1/2, 0, 1/2) and will be seen to arise 
from magnetic scattering. Two unusual aspects of the 
magnetic diffraction pattern will be shown to define the 
magnetic structure: (1) The observable magnetic peak 
with the shortest wave vector is at Q = 1.75 A~^ and 
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FIG. 1: Sr2F2Fe20S2 (a) neutron diffraction pattern for T = 
300 K (red - measured, black - calculated, green - difference 
between the measured and calculated patterns, blue vertical 
lines - calculated nuclear peak positions) and (b) the unit cell, 
with the red arrows indicating the Fe spins orientation in the 
AFM state, and (c) the ab plane view of the Fe2 0S2 layer. 



is associated with indexes (0.5, 0, 4.5) and (0.5, 1, - 
0.5). This points to a longitudinally modulated mag- 
netic structure where the polarization factor extinguishes 
peaks at wave vectors that do not have an adequate com- 
ponent transverse to the in-plane direction of unit cell 
doubling (here the a-direction) ; (2) The (0.5, 1, -0.5) 
magnetic Bragg peak stands out by a sharp leading edge 
which is not apparent at (0.5, 1, 1.5). This indicates 
strong correlation length anisotropy with much longer 
correlations in the a-b plane than along c. 

For a systematic analysis of possible magnetic struc- 
tures we used representation analysis and Rietveld re- 
finement as implemented in SARAiP^ and fullprof.'A^ 
For the I/Ammm symmetry with a single propagation 
vector of k = (1/2, 0, 1/2), the iron sites are divided 
into two independent sets: those forming Fe-O-Fe chains 
along the in-plane a-component of k = (1/2, 0, 1/2) and 
those forming Fe-S2-Fe chains along a (Fig. [ij;). Using 
Kovalev notation,'^^! the representation Tjnag associated 
with such magnetic structures decomposes into two irre- 
ducible representations (irreps), Fi and Fa in the form of 
Fmag = Fi -|- 2F3, with all three basis vectors projected. 
Given the observations above which are quantitatively 
born out by the Rietveld analysis, F3 with moments along 
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FIG. 2: Neutron diffraction patterns for T = (a) 300 K, (b) 
100 K and (c) 12 K and 20 K (averaged), showing the devel- 
opment of magnetic peaks below Tjv = 106(2) K (green tick 
marks). The magnetic peaks are refined with a Warren-like 
line shape (red lines, see text) that consistently accounts for 
the magnetic diffraction data. The fitting background and 
the difference between the measured and calculated patterns 
are shown in black dotted lines and brown solid lines, respec- 
tively. The blue dotted line represents the best fit profile 
assuming an isotropic magnetic correlation length, which is 
however inconsistent with the experimental data. 



the in-plane component of k is the only irrep that is con- 
sistent with the data. 

As the refinement based on the model of isotropic cor- 
relation length (blue dotted line. Fig. [2j:) failed to ac- 
count for the experimental data (x^ = 7.496), the correla- 
tion length anisotropy was quantified by replacing in the 
Rietveld analysis of magnetic diffraction the spherical av- 
erage of d^{Q — t) by the spherical average of the follow- 
ing normalized anisotropic correlation function (Warren- 
like function^*'): 
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The gaussian and lorentzian respectively describe long 
range in-plane correlations and short range correlations 
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FIG. 3: (a) DC magnetization measured for H = 1 T. The 
ZFC and FC data are shown as sohd and open symbols, re- 
spectively; inset: the hysteresis loop M{H) for T — 2 K. (b) 
The T-dependent staggered magnetization obtained by fitting 
diffraction data. The solid line shows the Onsager solution of 
the 2D square lattice Ising model, (c) semi-log plot of resis- 
tivity vs. 1/T (symbols), with a linear fit that gives Eg = 
0.28 eV. The raw p{T) data are shown in the inset; (d) zero 
field heat capacity. After subtracting the phonon contribution 
(dashed line), the magnetic heat capacity Cm{T) is plotted vs. 
t — \T — Tn\/Tn (inset), showing logarithmic divergence near 
Tn for both T < Tn (solid) and T > Tn (open). 



along c. The polarization factor associated with the 
structure was implemented within the spherical average. 
The corresponding fit shown in Fig. [2] provides an ex- 
cellent account of the magnetic diffraction data (x^ = 
1.033). The in-plane peak width cr± = 0.0036 1"^ is 
near the resolution width extracted from the nearby nu- 
clear peak: ares = 0.0023 A^^, indicating the in-plane 
correlation length exceeds 310 A. By contrast, the cor- 
relation length along c, — 17(3)A, is indistinguish- 
able from the lattice parameter c, so the low temperature 
magnetic order can rightfully be characterized as quasi- 
two-dimensional. These results are consistent with the 
initial phenomenological observations. 

A well-known ambiguity presents itself in interpreting 
the results. The Fe-O-Fe and Fe-S2-Fe chains associ- 
ated with ki — (1/2, 0, 1/2) have identical, real magnetic 
structure factors for allowed magnetic Bragg peaks in the 
coUinear Fa spin structure. The magnetic diffraction in- 
tensity from a single wave vector domain thus depends 
only on the sum of the ordered moment on the two sites. 
By calibrating the magnetic diffraction intensity against 
the nuclear diffraction intensity in the conventional fash- 
ion, the magnitude of the total moment comes out to 
be 4.7(1) ^B- It is not possible, on the basis of diffrac- 
tion alone, to distinguish the collinear structure from a 
non-collinear two-k structure where perpendicular Fe-O- 
Fe chains are modulated antiferromagnetically (Fig. [ij;) . 
However, while the collinear structure defines two dis- 




FIG. 4: (a) Non-magnetic DOS and (b) Fermi surface of 
Sr2F2Fe20S2, showing a narrow Fe 3d band and quasi-two- 
dimensional features. 



tinct magnetic sites with an average moment of 2.35(5) 
(=4.7(l)/2) /is/Fe^+, the two-k structure implies an or- 
dered moment of 3.3(1) (= 4.7(1)/V2) fis/Fe^^ on all 
sites. Of these, only the latter structure yields a mo- 
ment comparable to the spin-only localized moment of 
4 /is/Fe^^, and is consistent with the observation of a 
single site Fe^"*" hyperfine split Mossbauer spectrum.HSl 

The noncollinear spin structure can be described as 
longitudinally polarized antiferromagnetic order on per- 
pendicular Fe-O-Fe chains that intersect at oxygen sites. 
As for the "2322" compounds, the intra-layer exchange 
interactions between the nearest and next nearest neigh- 
boring Fe sites involve three distinct exchange paths 
(green arrows. Fig. [ij;): Ji links spins within chains 
via the 180° Fe-O-Fe super-exchange interaction and is 
expected to be dominant. J2 links perpendicular chains 
through the 90° Fe-O-Fe and the approximately 90° Fe- 
S-Fe super-exchange interactions. J3 corresponds to the 
Fe-S-Fe super-exchange with a bond angle of 100.22°, 
which links parallel spin chains. The observed magnetic 
structure is favored by AFM Ji and FM J3 interactions. 
Isotropic J2 interactions are frustrated due to the 90° 
angle between nearest neighbor spins. 

The proposed magnetic structure differs from those 
inferred for the structurally analogous "2322" systems 
LaaOsTzSez (T = Mn, Fe, Co), for which the T ^ Mn 
compound shows a G-type AFM structure with all spins 
parallel to the c axis,'^ the T = Fe compound exhibits a 
bi-coUinear AFM structure with all spins parallel to the 
a axis,'J^ while for T = Co a non-collinear plaquette AFM 
structure with 90° angles between the spin orientations 
of neighboring Co is observed.'^i' None of these structures 
is consistent with the magnetic diffraction pattern for 
Sr2F2Fe20S2 (Fig. 2c). The unusual diversity of mag- 
netic structures indicates frustration in the magnetism of 
the checkerboard Fe2 0X2 spin lattice. 

The lamellar magnetism revealed by the neutron 
data is also reflected in the thermodynamic properties. 
The DC magnetization M(T), measured in a Quantum 
Design Magnetic Property Measurement System (QD 
MPMS) for = 1 T, shows a peak at Tjv = 106(2) 
K (Fig. [3^), as expected for long range AFM order- 
ing. While the Neel temperature is consistent with a 
previous report ,1^^ the anomaly reported here is much 
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TABLE I: Relative energies AE (meV per unit cell) of six spin configurations and exchange constants Ji (meV/bond) for 
different values of U . 
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sharper, indicating, as is generally the case in frustrated 
magnets, that sample quality can affect the transition. 
Heat capacity data (Fig. [Sji), measured on a QD Phys- 
ical Property Measurement System (PPMS) also shows 
a sharp anomaly. By subtracting the estimated phonon 
contribution from a polynomial fit of the Cp{T) data at 
temperatures away from T^r (dashed line. Fig. |3ji), the 
singular magnetic specific heat Cm was extracted and 
plotted vs. t = |r - Tn\/Tn in the inset of Fig. |3}i. 
The observed logarithmic divergence of Cm is expected 
for a two dimensional Ising system. Fig. [Sfj shows the 
temperature dependence of the sublattice magnetization 
m(T) inferred from diffraction. Consistent with the spe- 
cific heat, the data resemble the Onsager result for the 
square lattice Ising model (solid line). 021 Ising criticality 
is consistent with strongly anisotropic exchange interac- 
tions which also characterize the iron pnictides.'^ 

The zero-field-cooled (ZFC) and field-cooled (FC) 
magnetization curves (full and open symbols, respec- 
tively, in Fig. [3^) exhibit irreversibility below Tjv and 
this is also reflected by the magnetic hysteresis loop for 
T = 2 K (inset, Fig. [3^). This indicates glassy be- 
havior, consistent with the short range inter-plane cor- 
relations. While commensurate AFM order on a square 
lattice comes in two time-reversed versions, the spins in 
each of the perpendicular sets of chains can be separately 
reversed in the two-k structure, leading to a total of 
four different versions. Noting substantial in-plane cor- 
relations above T/v indicated by the susceptibility max- 
imum around 150 K (Fig. [3|, consider now the stacking 
of such spin planes along c. To establish long range or- 
der along c, Dzyaloshinskii-Moriya (DM) interactions be- 
tween spins displaced by c/2 and exchange interactions 
between spins displaced by c, which - depending on the 
sign of the DM interaction - could even be frustrated, 
must select just one out of these four copies. Consider- 
ing the slow dynamics associated with transitioning be- 
tween the four different Ising spin planes, it is plausible 
the process fails due to kinetics and/or disorder. 

Like aU other "2322" and "22212" compounds, 
Sr2F2Fe2 0S2 is an insulator for which the temperature 
dependent resistivity p{T) (inset. Fig. |3j:) can be de- 
scribed by p oc e^a/^sT ^ shown in the Arrhenius plot 
(symbols. Fig. |3j:). A linear fit of Inp vs. \/T (solid line, 
Fig. [3b) suggests an activation gap Eg — 0.28 eV, which is 
considerably larger than the previously reported value of 
0.10 eV but equal to Eg = 0.28 eV for Ba2F2Fe20Se2.'i2l 
At temperatures lower than 250 K, the resistivity starts 
to deviate from the activation gap model, which suggests 
the importance of hopping between localized gap states, 



as observed in LaMnP0i_j,Fa;.'2il 

To connect our experimental observations with the 
electronic structure of Sr2F2Fe20S2, DFT calculations 
were performed using the full-potential linearized aug- 
mented plane wave (FP-LAPW) method implemented in 
the wien2k package The generalized gradient approx- 
imation (GGA) was used for the exchange correlation po- 
tential, and different effective onsite Coulomb repulsion 
values U from 1.5 to 4.5 eV were considered. An Fe 2>d 
band between -2 eV and 1.5 eV is revealed in the non- 
magnetic density of states (DOS) spectrum (Fig. [4k). 
The 3d band is narrower than that of LaOFeAa^S and 
BaFe2As2,'^a direct consequence of the larger Fe square 
lattice in Sr2F2Fe20S2 and the quasi-one-dimensional 
connectivity. Since the band width is proportional to the 
kinetic energy the narrow id band suggests enhanced 
electron correlations U/t, leading to the insulating be- 
havior and a relatively high Neel temperature. The Sr-F 
layers (red and dark yellow lines. Fig. |4^) contribute 
negligible DOS to the Fermi level, hence the inter-layer 
electron hopping is strongly reduced, giving rise to the 
quasi-two-dimensional electronic structure. The low di- 
mensionality of Sr2F2Fe20S2 is also reflected in the Fermi 
surface, which mainly consists of quasi-2D ribbons ex- 
tending along the F — Z direction (Fig. ^p) . 

The magnetic ground state and exchange constants Ji , 
J2 and J3 are estimated from spin-polarized calculations 
for the following six spin configurations (as shown in Ref. 
2): (1) FM, (2) checkerboard AFM (AFMl), (3) single 
stripe-l-FM (AFM2), (4) double stripe-|-FM (AFM3), (5) 
(1/2, 1/2) state with antiparallel alignment across the 
O sites (AFM4), (6) stripe AFM (AFM5). For differ- 
ent values of the Coulomb repulsion U , the ground state 
energies are listed in Tabic |T| 

For ?7 = 0, finite DOS exists at the Fermi level for all 
spin configurations, while a Mott gap opens for finite U . 
Out of all the calculated spin configurations (Table [l| , 
AFM4 is found to be the lowest energy state, with the 
ordered moment determined to be 3.49 /is/Fe^"*". Al- 
though the specific orthogonal spin configuration is not 
uniquely identifiable in the DFT results due to software 
limitations, the moment value is in very good agreement 
with the value associated with the two-k magnetic struc- 
ture (Fig. [l|;). In addition, the signs of the exchange 
coupling constants Ji and J2 (Table |l]) suggest an AFM- 
type Fe-O-Fe superexchange coupling and an FM-type 
coupling along the Fe-S-Fe paths, which are consistent 
with the two-k structure rather than the collinear struc- 
ture with ki ~ (1/2, 0, 1/2). As suggested by the phase 
diagram proposed in Ref. [H the condition for AFM4 
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to be the ground state is Ji > J2 and J3 < 0, which 
also agrees well with the calculated exchange constants 
of Sr2F2Fe20S2. By calculating the energy difference 
between parallel and antiparallel spin configurations, the 
inter-layer interaction J± is estimated to be three orders 
of magnitude weaker. 

In conclusion, we have provided evidence for a non- 
coUinear magnetic structure in the layered oxychalco- 
genide Sr2F2Fe2 0S2. The Fe spin lattice consists of 
orthogonally intersecting ID AFM chains, with frus- 
trated near neighbor interactions - a structure that to our 
knowledge has not previously been observed in a square 
lattice Mott insulator. The phase transition is of the 
2D Ising variety and electronic two-dimensionality, inter- 
layer magnetic frustration, and four degenerate versions 
of the two-k magnetic structure are contributing factors 



to magnetic hysteresis and short-range interlayer corre- 
lations. The intersecting chain-like structure is the basis 
for the unusual non-collinear magnetism in Sr2F2Fe20S2 
and through effective electronic one-dimensionality per- 
haps also its Mott insulating state. The variety of mag- 
netic structures suggests significant potential for tuning 
electronic correlations in the extended "2322" family of 
materials. 
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